&F  TECHNICAL  REPORT  NO.  5790 


July  1949 


TRANSONIC  FLOW 

THROUGH  NEARLY  CIRCULAR  DUCTS 
WITH  GRADUALLY  VARYING  CROSS  SECTION. 

Don  Blose 


distribution  statement  a 

Approved  for  Public  Release 
Distribution  Unlimited 


UNITED  STATES  AIR  FORCE 
AIR  MATERIEL  COMMAND 

W RIGHT -PATTE RSON  AIR  FORCE  BASE,  DAYTON,  OHIO 


Reproduced  From 
Best  Available  Copy 


20000515  045 


NOTE 

When  drawings,  specifications,  and  other  data  prepared  by  the 
War  Department  are  furnished  to  manufacturers  and  others  for  use 
in  the  manufacture  or  purchase  of  supplies,  or  for  any  other  pur¬ 
pose,  the  Government  assumes  no  responsibility  nor  obligation 
whatever;  and  the  furnishing  of  said  data  by  the  War  Department 
is  not  to  be  regarded  by  implication  or  otherwise,  or  in  any  manner 
licensing  the  holder,  or  conveying  any  rights  or  permission  to 
manufacture,  use,  or  sell  any  patented  inventions  that  may  in  any 
way  be  related  thereto.  ' 

The  information  furnished  herewith  is  made  available  for  study 
upon  the  understanding  that  the  Government's  proprietary  interests 
in  and  relating  thereto  shall  not  be  impaired.  It  is  desired  that  the 
Patent  &  Royalties  Section,  Office  of  the  Judge  Advocate,  Air 
Materiel  Command,  Wright-Patterson  AFB,  Dayton,  Ohio,  be 
promptly  notified  of  any  apparent  conflict  between  the  Govern¬ 
ment's  proprietary  interests  and  those  of 


-  Espionage  Act '  *  i-  v 

Notice:  “This  document  contains  information  affecting  the 
national  defense  of  the  United  States  within  the  meaning  of  the 
Espionage  Act,  50  U.  S.  C.,  31  and  32,  as  amended.  Its  transmission 
or  the  revelation  of  its  contents  in  any  manner  to  an  unauthorized 
person  is  prohibited  by  law.” 

“The  above  Espionage  Notice  can  be  disregarded ^unless  this 
document  is  plainly  marked  with  a  security  classification  as 
“Restricted,”  “Confidential,”  “Secret,”  or  “Top  Secret.”^.;:  .  ; 


The  U.  S.  Government  is  absolved  from  any  litigation  which 
may  ensue  from  the  contractor’s  infringing  on  the  foreign  patent 
rights  which  may  be  involved. 


REPORT  DOCUMENTATION  PAGE 


Form  Approved 
OMB  No.  0704-0188 


1.  AGENCY  USE  ONLY  (Leave blank) 


the  collection  of  information,  send  comments  renaming  mis  uuraen  estimate  w  any  umei  m  . . . - ^  — n  D.  moot  nr  onsm 

Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington,  VA  222024302,  and  to  the  Office  of  Management  and  Budget,  Paperwork  Reduction  Project  (0704-0188),  Washington,  DC  20503. 

2.  REPORT  DATE  3.  REPORT  TYPE  AND  DATES  COVERED 

JULY  1949  FINAL  JANUARY  -  JULY  1949 

OTTLE  AND  SUBTITLE - 5.  FUNDING  NUMBERS 

TRANSONIC  FLOW  THROUGH  NEARLY  CIRCULAR  DUCTS  WITH 
GRADUALLY  VARYING  CROSS-SECTION 

6.  AUTHOR(S)  ' 

DON  BLOSE 


7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

AIRCRAFT  LABORATORY 

ENGINEERING  DIVISION  (WIND  TUNNEL  BRANCH) 
AIR  FORCE  MATERIEL  COMMAND 
WRIGHT-PATTERSON  AFB,  OH  45433 


9.  SPONSORINGfMONITORING  AGENCY  NAMEIS)  AND  ADDRESS(ES) 

AIRCRAFT  LABORATORY 

ENGINEERING  DIVISION  (WIND  TUNNEL  BRANCH) 
AIR  FORCE  MATERIEL  COMMAND 
WRIGHT-PATTERSON  AFB,  OH  45433 


ITT.  SUPPLEMENTARY  NOTES 


8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 


10.  SPONSORINGfMONITORING 
AGENCY  REPORT  NUMBER 

AF-TR-5790 


12a.  DISTRIBUTION  AVAILABILITY  STATEMENT 

APPROVED  FOR  PUBLIC  RELEASE;  DISTRIBUTION  UNLIMITED 


12b.  DISTRIBUTION  CODE 


13.  ABSTRACT  (Maximum  200  words)  _  TT,™ 

A  METHOD  FOR  THE  ANALYSIS  OF  A  TRANSONIC  FLOW  FIELD  IN  A  NEARLY  CIRCULAR  DUCT  WITH 
GRADUAL  CHANGES  IN  CROSS-SECTION  IS  DEVELOPED.  A  LINEARIZED  DIFFERENTIAL  EQUATIONF  FOR 
THE  DEVIATION  FROM  THE  RESULTS  OF  THE  ONE-DIMENSIONAL  THEORY  IS  PRESENTED  AND  IS  SOLVED 
BY  THE  SUPERPOSITION  OF  PARTICULAR  SOLUTIONS  OBTAINED  BY  A  PRODUCT  HYPOTHESIS.  THE 
POTENTIAL  EQUATION,  IN  THIS  CASE,  WAS  SIMPLIFIED  FOR  THE  VICINITY  OF  THE  SONIC  VELOCITY. 
APPLICATION  TO  A  CIRCULAR  SUBSONIC-SUPERSONIC  DUCT  IS  MADE  IN  ORDER  TO  DETERMINE  THE 
MAGNITUDE  OF  PROPAGATED  SUBSONIC  DISTURBANCES  AND  THE  RESULTING  VELOCITY  DISTRIBUTION. 


14.  SUBJECT  TERMS 

TRANSONIC  FLOW,  SUBSONIC  DUCT,  SUPERSONIC  DUCT 


15.  NUMBER  OF  PAGES 

36 

16.  PRICE  CODE 


17.  SECURITY  CLASSIFICATION 
OF  REPORT 

UNCLASSIFIED 


18.  SECURITY  CLASSIFICATION 
OF  THIS  PAGE 

UNCLASSIFIED 


19.  SECURITY  CLASSIFICATION 
OF  ABSTRACT 


UNCLASSIFIED 


20.  LIMITATION  OF 
ABSTRACT 


Standard  Form  298  (Rev.  2-89)  (EG) 

Prescribed  by  ANSI  Std.  239.18 

Designed  using  Perform  Pro,  WHS/DIOR,  Oct  94 


AF  TECHNICAL  REPORT  NO.  5790 


TRANSONIC  FLOYf 

THROUGH  NEARLY  CIRCULAR  DUCTS 
WITH  GRADUALLY  VARYING  CROSS-SECTION 


Don  Blose 


July  191+9 


Published  by 
UNITED  STATES  AIR  FORCE 
AIR  MATERIEL  COMMAND 


Wri ght-Pat t e rs on  Air  Force  Base,  Dayton,  Ohio 


WF-(A).0-21  OCT  49 


42 


ABSTRACT 


A  method  for  the  analysis  of  a  transonic  flow  field  in 
a  nearly  circular  duct  with  gradual  changes  in  cross-section 
is  developed.  A  linearized  differential  equation  for  the 
deviation  from  the  results  of  the  one-dimensional  theory  is 
presented  and  is  solved  by  the  superposition  of  particular 
solutions  obtained  by  a  product  hypothesis.  The  potential 
equation,  in  this  case,  was  simplified  for  the  vicinity  of 
the  sonic  velocity. 

Application  to  a  circular  subsonic-supersonic  duct  is 
made  in  order  to  determine  the  magnitude  of  propagated  sub¬ 
sonic  disturbances  and  the  resulting  velocity  distribution. 
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SECTION  I 


INTRODUCTION 

This  report  presents  an  analytical  method  for  the  calculation  of 
the  flow  field  in  a  nearly  circular  duct  with  relatively  gradual  changes 
in  cross-section.  Due  to  simplifications  introduced  by  means  of  the 
transonic  law  of  similarity*  the  method  in  its  present  form  would  be 
restricted  to  flow  fields  with  Mach  numbers  close  to  unity. 

The  investigation  was  initiated  in  connection  with  a  plan  for 
modifying  the  10-Pbot  Wind  Tunnel  at  Wright-Patterson  Air  Force  Base. 

It  was  intended  to  vary  the  cross-sectional  area  with  small  inserts* 
rather  than  by  a  variation  of  the  radius.  A  method  was  evolved  to 
obtain  the  magnitude  of  deviations  from  the  ideal  parallel  flow  which 
might  be  caused  by  such  an  insert  design.  Before  the  actual  computa¬ 
tions  were  made  the  design  was  modified.  However*  the  method  itself 
is  of  interest,  in  that  there  are  many  other  applications.  An  appli¬ 
cation*  referring  mainly  to  a  modified  test  section  design*  will  be 
given  in  the  second  part  of  this  report. 

The  present  analysis  begins  with  a  basic  flow  as  given  by  one- 
dimensional  theory.  Any  deviations  from  the  basic  flow  are  considered 
small,  so  that  they  can  be  determined  by  a  linearised  partial  differen¬ 
tial  equation.  The  general  solution  to  the  differential  equation  is 
expressed  by  the  superposition  of  particular  solutions*  each  of  which 
is  the  product  of  a  Bessel  function  dependent  on  the  radial  distance 
y  from  the  duct  axis ,  a  trigonometric  function  of  the  angular  position 
<o*  and  a  function  of  the  axial  position*  x.  Of  special  interest  is 
the  last  function  which  is  influenced  by  the  duct  shape. 

If  one  is  interested  only  in  supersonic  flow  analysis,  an  alternate 
method  of  computation  is  the  well-known  method  of  characteristics  •  The 
method  Of  characteristics  is  more  general*  in  that  the  condition  of  a 
slow  change  in  cross-section  is  not  required.  The  method  of  character¬ 
istics  for  three  dimensions  is  similar  to  the  present  method  in  the 
requirement  of  a  nearly  circular  cross-section.  The  major  advantages 
of  the  present  method  are  that  the  subsonic  influences  can  be  deter¬ 
mined  and  that  being  essentially  more  analytical*  a  better  picture  of 
the  flow  behavior  in  general  cases  may  be  obtained.  In  some  cases* 
depending  upon  the  duct  shape  *  the  solutions  can  be  obtained  in  a  - 
less  tedious  manner. 

Briefly*  the  report  starts  with  the  simplification  of  the  general 
flow  differential  equation  by  means  of  the  transonic  law  of  similarity. 
Then  a  general  method  of  solution  is  developed*  and  finally  this  result 
is  applied  to  some  specific  example . 
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SECTION  II 


THEORETICAL  CONSIDERATIONS 

Simplification  of  the  Plow  Differential  Equation  Near  the  Sonic  Speed 

This  analysis  is  based  on  a  cylindrical  coordinate  system  (Figure  1), 
Let  the  x  axis  be  the  duct  center  line,  y,  the  radial  distance  from  the 
x  axis,  and  <w  ,  the  angular  position  measured  clockwise  frcsn  the  vertical 
position.  The  velocity  potential  is  denoted  by  $  .  Partial  derivatives 

are  denoted  by  subscripts  (i.e.,  <J>^) 

The  general  three-dimensional  equation  for  a  steady  potential  flow 
(derived  in  Appendix  A)  is  given  by 

U.U  ti  t  4^.0 

If  deviations  from  a  parallel  sonic  flow  are  small,  the  simplified 
transonic  equation  results.  (Derivation  may  be  found  in  Appendix  B).  • 

(1.2)  -(<+>)  <P»«  q>« -«-<?„+  +  =  0 

where  ^  -  S, 

r* 

y  =X 

U-3) 

and  rQ  is  the  radius  at  the  throat. 

The  Pressure  and  Density  Relations 

The  new  potential  term  <P(x,y)  is  directly  related  to  the  velocity. 

However,  one  may  be  interested  in  the  pressure  and  density  relations  at 
various  positions  along  the  duct.  Bernoulli's  law  for  compressible  non- 
viscous  gases  is  given  by 

(1.1*)  ^  tU4UtVAV:0 

The  last  term  is  negligible,  based  on  earlier  assumptions . 

Because  Alls  C?*  f  one  can  write 

*  * 

Since  ,  the  preceding  equation  is  rewritten  as 
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0-5)  ^  Cpx 

The  Boundary  Conditions 

Let  7  be  the  value  of  y  at  a  point  on  the  duct  surface  so  that 
r  ■  r  (x,  o>).  This  function,  in  the  case  of  a  vertical  line  of  symmetry 
is  expressed  as  ^ 

(1.6)  r  =  r,  U)  -h  2  (f*  (x)  cos  «ui) 


The  functions  r0,  r^,  r2, . rn  are  obtained  by  a  Fourier  analysis 

at  the  specific  x  positions,  and  are  written  as 


— 

(1.6a)  V0(X)=^  J  <i 


CO 


(i.6b) 


cos  n co  du» 


Let  r  ■  — —  . 

r0 

The  boundary  condition  of  aero  flow  normal  to  the  duct  surface  must 
be  fulfilled.  If  the  differences  from  the  sonic  velocity  are  small, 
this  condition  is  expressed  as 

(1.7)  -  ^*<PY=© 

(The  derivation  may  be  found  in  Appendix  C). 

Let  us  introduce  for  the  slope 

&•»> 

Then  one  has  the  relation 

(1.9)  <Py  =  I  «-,«•*  V«» 

Where  "©*v  is  an  additional  slope  for  the  unsymmetrical  configurations . 

The  Basic  Hypothesis  for  the  Solution 

The  solution  is  based  on  one-dimensional  flow  theory,  in  which  the 
velocity  distribution  over  each  cross  section  is  considered  constant. 

The  cross  section  is  written  as  functions  of  the  distance  from  the 
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minimum  cross-section.  The  potential  of  this  basic  flow  is  denoted  by 

,  where  —  QPe>(  =.  -p(r,r0V  Specifically,  as  derived  in  Appendix 

D,  this  gives 

(1.10)  <PB<  =  (£•,**)* 

The  entire  potential  is  expressed  by 

(1.11)  Cpr  Cpe(y:H 

The  CP  represents  the  deviations  from  the  basic  flow,  and  shall  be 
considered  small  enough  that  higher  order  terms  can  be  neglected.  If  the 
values  of  the  complete  potential  (Equation  (1.11))  is  inserted  into  the 
transonic  differential  equation  (Equation  (1.2))  and  second  order  terms 
of  Cp  are  neglected,  the  following  inhomogeneous  equation  is  given* 

(1.12)  V W*”) <*». ‘P*.. 

For  the  representation  of  ,  let  us  introduce  the  following  functions. 
Let  a^  Bessel  function  of  V  order  and  of  the  argument  1]  be  denoted  try  Jv ('»)). 
Let  vjv).  be  the  ^ argument  (where  y.  -1,2,3  .  .  .)of  the  Bessel 

function  Jv  ,  for  which  4  2x1a),  o  .  The  functions  2?Vyk,  are  given  by 

r 

(1.13)  ZV(4<y)=  J.fy  yj,r) 

The  functions  ja.  fulfill  the  following  differential  equation 

d.ii.)  (Z,r)±  &&£>  +  C-yjV-  v?  )  i  *>.  -  O 

There  exists  a  relationship  of  orthogonality  between  the  2VUi  ,  so 
that  for  v  ' 

(1.15)  j  Y  iyy.  dy-O 
Let  us  introduce 

(1.16)  (  Y  <*y-Kv>* 

Jo 

A  hypothesis  for  the  perturbation  velocity  potential  is  chosen  in  such  a  way 
that  the  boundary  conditions  are  fulfilled* 
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(1.17) 


<p  -  t  v1  e.(*)  +  5  (i  Y*  ^  cos 

+  I  (a«>W  +  |,  cos,,u) 

/*r\  r 

The  functions  and  aVj».(*)  (  V  ■  1*  2,  3  •  •  .  )  must 

be  determined  in  such  a  way  that  the  differential  equations  for  q> 

(Equation  (1.12)  is  fulfilled.  Inserting  this  hypothesis  into  Equation 
(1.12),  one  obtains 

(1.18)  +  2  0i  £vr  cos  <Pb« 

I,  I,  *»»*  c“  v“>3  *j  I,  cv“ vu 

+  -  7.  Zv^. eo* vu)  -('*')[  i,  .  v” 

cos  <“(<?.,  ■©•»  *  <pgt,  e>j] = 


ail  the  terms  of  this  equation  which  contain  trigonometric  terms  of 
(where  V  -  1,  2,  3,  h,  etc.)  must  vanish  separately. 

For  V  -  0  one  has  the  condition 

-<<♦»)[$«*  (h2b^+  I 

(1.19)  <X.^  (  f  -  ft*')  C  + 

—  -2-©^-  0$*0  ^B*  CpBun 


are 


The  conditions  for  V  ■  1,  2,  3>  •  •  ^ 

-OH-Of/Je  f  Z,,.V<Po„  J  (o-^ivY) 

*  Z4*1 

*!,  ^  - 1 4r  -  £ z*^ 

"  (3r)  V*  (^By^v  +  ^B**  ^v) 

It  can  readily  be  seen  from  Equation  (1*1.0)  that 

(1.21)  e«.:  -(*+>)  q>B»  <Pb». 
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Further  simplifications  are  made  using  Equation  (1.11+)  with  Equa¬ 
tions  (1.19)  and  (1.20),  respectively. 

(1.22) 

“1,  KV  vl  z  (2r  (*•*  <  -  $B**  < ) 

(1.25)  L^v/*(^Pb»  av^  ■'*  ^p©**  avjOj 

-  J-,  IX  ^  av^]  =  7V  C<Pox <pBu¥  -ei) 

The  equations  for  the  perturbation  functions  cx.  v^».  are  obtained 
by  multiplying  with  the  term  y  2rvj*.  and  integrating  between  y  ■  0  to 
y  ■  1.  Because  of  the  orthogonality  relations  for  the  Bessel  functions, 
all  with  the  terms  J*  y  2. y/<  Jy  drop  out  except  for  . 

Inserting  these  relations  in  Equations  (1.22)  and  (1.23)  one  obtains 
Equations  (1.21+)  and  (1.25),  respectively. 

U.2W  Vt**  °v *  cpa»«  +  %- 

~  H  .j*  "t  +  ”®0 

where 

(1.25)  K^-  fs.,w  y3 -Ay  =  h  1.(7  V) 

(1.26)  iv=  f  s.„ (7)  7  * y  =  T  V<v 

^  O 

(Derivation  may  be  found  in  the  Appendix  E).  Also, 

(1-27)  l  q>8,  a,;  +  cpB,„  av;  *  5k  avr] 

"  “  Kij»  "T  C  "*■  ^Bn*  ^  ") 

(1.28)  X,j»  3  f  ’  £.,.(?)  7  <ly  =  X(:V  7> 

(1.29)  K,^=££vr(7>)  7"'  4y-^.^v(v  +") 

From  Equations  (1.21+)  and  (1.27),  and  the  boundary  conditions  which 
will  be  given  later,  the  values  of  o-^can  be  found.  The  resulting 
expression  for  the  perturbation  potential  (Equation  (1.17))  is  written  in 
a  more  convenient  manner  if  one  expresses  the  l/2  y2  and  x  yV  by  means 
of  Z^and  2vr.  .  ( 


where 


<=£  K± 

2-  X 


Let  us  introduce 


°r- 


(1.30) 

(1.31) 


bv^  +  ^  v ~  V  )> 


The  perturbation  potential  is  then  based  entirely  in  terns  of  the 
functions 


2.  v  so  that 

(1.52)  v  «  £  C  <4  2-rl  *  J.  I,  C  a»r2T] 

*- 

One  can  then  write  the  differential  equations  for  ty  the 

introduction  of  Equations  (1.30)  and  (1.31)  into  Equations  (1.24)  and 
(1-27)- 

_  *  //  _  *  /  *  K®>» 

^B*  °"*>*  +  <Pb**  aor+  'JS  +  l  •>*  “ 

*  "  -  -♦  '  -  ik  -  Ki 


(1.33) 


(1.54)  +  <$>©*,  CXv/k  +  <XV^  =  ^  -©v 


CP  and  its  first  derivatives  with  respect  to  x  are  continuous 
Therefore,  according  to  thi3  relation,  ~  J  J,x — '  “  J — 5 
are  continuous. 


Ow~  ^  and  their  first  derivatives 


One  might  ask  why  the  hypothesis  (Equation  (1.32))  was  not  used 
from  the  beginning  instead  of  Equation  (1.17).  The  reason  is  that  the 
second  derivative  with  respect  to  y  does  not  converge  (because  of  the 
character  of  the  function  yv  )•  Therefore,  all  of  the  considerations 
regarding  the  fulfillment  of  the  differential  equation  for  would 
have  been  meaningless. 


Method  of  Solution 

♦ 

The  solution  of  the  perturbation  functions  for  the  differen¬ 

tial  equations  (1.33)  and  (1.34)  is  found  for  ary  duct,  within  the  defined 
limitations,  by  means  of  numerical  integration.  For  suitable  configurations, 
analytical  solutions  can  be  found.  In  some  cases,  a  piecewise  analytical 
solution  may  be  satisfactoiy.  In  any  case,  the  pressure  distribution 
along  the  axis  is  influenced  only  by  the  expressions  with  V  ■  0,  since 
the  nature  of  higher  order  Bessel  functions  is  such  that  there  is  no 
effect  along  the  axis. 

Deviation  from  a  parallel  flow  are  found  from  <P^  *  The  only  expres¬ 
sions  which  are  of  Importance  along  the  axis  are  those  for  V  ■  0  and 
V  ■  2.  Actually,  an  infinite  number  of  terms  are  needed  to  describe 
the  flow  completely,  but  the  higher  terms  in  y-  and  y  are  very  small 
and  are  of  influence  only  close  to  the  wall. 
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Conditions  for  the  Solution  of  Equations  (1.33)  and  (1.34) 


Equations  (1.33)  and  (1.34)  are  of  the  second  order,  and  consequently 
the  solutions  must  be  determined  by  two  boundary  conditions.  The  first 
condition  is  obvious.  With  the  assumption  of  a  subsonic  duct  extending 
to  negative  infinity,  one  may  assume  that  the  values  of  the  perturbation 
terms  <xv^  should  vanish  as  one  approaches  negative  infinity. 

The  character  of  the  second  condition  is  less  obvious.  One  may 
first  ask  at  which  point  such  a  condition  might  be  given.  Certainly  it 
cannot  be  a  point  in  the  supersonic  region,  for  the  fulfillment  of  the 
condition  will  affect  the  entire  solution,  including  points  upstream  of 
the  point  where  the  boundary  condition  is  given.  Since  this  would 
violate  the  law  of  the  forbidden  signals,  the  location  of  the  other 
condition  should  be  expected  to  be  in  the  subsonic  region.  Actually, 
mathematical  reasoning  leads  to  the  location  of  the  boundary  condition 
at  the  narrowest  cross-section.  If  one  rewrites  Equations  (1.33)  and 

(1.34)  in  such  a  way  that  the  coefficient  of  the  highest  derivative  of 
the  term  is  one,  the  coefficients  of  the  other  terms  become  singular 
at  the  sonic  velocity.  The  solution  of  the  differential  equation  may  be 
singular  at  any  point  where  these  coefficients  are  singular. 

Consider  now  only  the  solutions  of  the  homogeneous  part  of  the  dif¬ 
ferential  equation,  since  the  influence  of  change  in  the  arbitrary 
constants  and  consequently,  the  influence  of  the  boundary  condition  Is 
limited _to  that  part.  Let  the  basic  flow  velocity  be  represented  by 

Then  the  two  linearly  independent  solutions  of  the  homogeneous 
differential  equation  are  represented  by 

(1.35)  =•  P  (*) 

(1.36)  oSr  =  p(*) 

where  P  denotes  a  power  series  of  x,  which  starts  with  the  absolute  term. 

As  the  term  approaches  zero  (at  the  throat)  we  have  «  >  °  .  Then 
the  expression  o»„r  ,  from  Equation  (1.36)  will  tend  to  infinity  as  X-**0  . 
Since  this  is  physically  impossible,  the  second  condition  is  that 
remains  finite  at  the  throat. 

To  summarize  the  boundary  conditions ,  the  flow  deviation  expressions 
vj*.  must  be  finite  at  x  ■  0  and  zero  as  i  + 

SECTION  in 

APPLICATIONS 

Flow  Through  A  Circular  Subsonic— Supersonic  Duct 

We  apply  this  theory  to  an  axially-symmetric  transonic  duct  configura¬ 
tion.  This  duct  is  similar  to  the  modified  10-Foot  Wind  Tunnel  test  section 


S 
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at  Viright-Patterson  ^r  Forca^Baa^^Tha 

(Section  A*  Figures  U  an  5)  duct  expands  conically  (Section  C). 

profile  (Section  B  ,  and  finaUy. lubsoiic  Junction 

The  configuration  is  g^ometri^y  dete^  W  angle  of  expansion  - 

of  the  x-vulocity  deviations 
fU  a  parallel  flow  through  the  duct  channel. 

to  the  differ®^  ^tg;}5^°(wr,*#.ichhoccl^e’tafthe9 
SaSST^-  *£  Lv  .  -d’furthe^re.  in  the  function 
which  occurs  to  the  righthand  term. 

to  Section  A,  the  duct  is  cylindrical  and  consequently  the  radius, 
slope,  and  terms  are 

(2.1)  r=rv-!e-  *£&■ 

(the  value  of  “n  b»  ““  from  th*  sub3*’uent  f0rOula  (2J,))' 


(2.2)  ’O’o-  O 

,  t  /  u  (for  section  A) 

(2.3)  <?,,*=(£,  c) 

■i<.  nr-nfila  of  Section  B  is  an  approximation  for  a  circular  arc 

s?ii«Sr= jjyrssrs 

byt 

(2.4)  r  + 


(2.5) 

(2.6)  <Pe*=*(&T*y 


(for  section  B) 


Section  C  starts  tangential  to  the  parabolic  Section  B  and  expands 
conically  at  the  angle  <*  .  The  shape  terms  are 


(2.7)  r  -  i  +  (**  -  ^ 

(2.8)  6*0  =  Of 

(2.9) 


(for  section  C) 


These  values  inserted  into  the  differential  equation  (1.28)  lead 
to  the  specific  differential  equation  for  each  section. 

Section  A 


*  "  _*■  —  o 

(2.10)  CpB*  aV*  +  ^+1  V*  ” 
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Section  B 

■4fr  "  I  -#■ 

(2.12)  *  CXor  +  CXDr  +  E,r  ^=Elr  X 

The  constants  are*  _  x 

TLa*. 

(2.12)  -  /5_  £s\t 

W+l  R  ) 

»W 

(2.13)  E  -  ^ 

/i.  fjNt 

(?♦»  R  j 

Section  C 

The  solutions  of  the  differential  equations  (2.10),  (2.11),  and 
(2.110  will  now  be  determined.  The  complete  solution,  in  each  case, 
is  the  sum  of  the  general  solution  of  the  homogeneous  portion  and  a 
particular  solution  of  the  non-homogeneous  equation.  The  constants 
will  be  determined  later. 

Section  A 


The  solution  of  Equation  (2.10)  is 


V* 


21 i — \l  * 

v  f0-  /v*iW 


(2.15)  <£  =c,K  ev  +  clu  e 


Section  B 


The  particular  solution  of  Equation  (2.11)  is 

(2.14)  a!,*  -  E<rx 


where 


(2.17)  ~ 


—  _ 


*_  Xm-V- 


2  R  !•» 


3sk 

8+1 


The  homogeneous  part  of  Equation  (2.11),  through  a  simple  trans¬ 
formation,  is  of  the  form 
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the  solution'*'  of  which  is 

(2.20) 


The  term  denotes  any  linear  combination  of  Bessel  functions  of 
order  p.  The  general  solution  of  the  homogeneous  part  of  Equation  (2.11) 
then  is  such  that  the  complete  solution  becomes 


(2.21) 


*r 


=  E  E.,.*  +  c,r  3.  |>  *4J  ♦  C.jXL2  Er  X’J 


Section  C 

A  particular  solution  of  Equation  (2.1i*)  in  Section  C  is  given  by 


(2.22)  Ck*  = 

In  order  to  find  the  general  solution  of  the  homogeneous-  part  of 

Equation  (1.28),  the  transformations  u  =  (* “  and  <^(>0=cl^ 

are  carried  out.  The  differential  equation  is  then  written 

(2-23)  *U<-^)+ru(£)  +  E«/*<1=0 

where 

(2.21*) 


L»te 


Equation  (2.23)  has  the  form 

(2.25)  y  +  b%y  =  o 

where  the  solution2  is  given  as 

(2.26)  y  -  £x  [t  **] 

Consequently,  the  solution  of  Equation  (2.23)  is 

(2.27)  °1  "  ^ 


^  Jahnke,  E.  and  Emde,  F.,  Tables  of  functions.  Fourth  Edition.  Dover 
Publications,  19l*5>  P»  11*7  • 

2  Ibid.,  p.  11*7. 
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If  one  substitutes  for  u  and  g(u),  and  writes  the  constants  for 
the  explicit  Bessel  function,  one  finally  obtains 

(2.26)  0*0^=  <*  +  )4  CO+  C^T-^Ce,)] 

(2.29)  E5^  {C*>.l.|Ce,)'Cc>.  T^(e,)] 

where 

(2.30)  e,:| 


Determination  of  the  Constants 

The  constants  C,^,  Ci^  ,  ,  ...  .  and  c<r  which  occur  in 

the  solution  for  the  functions  o.«*  ,  are  found  by  equating  the  solutions 
of  a.**  ,  given  by  Equations  (2.15),  (2.21),  and  (2.26)  and  their  first 
derivatives  for  the  points  of  junction.  Furthermore,  there  are  the 
conditions  that  remain  finite  at  the  throat  and  that  they  -vanish 

as  x  tends  to  -  oo  . 

4fp 

The  condition  that  tends  to  zero  for  ,  according  to 

Equation  (2.15)  yields  Immediately 


(2.31)  C^-O 

The  condition  that  o f*.  is  finite  at  the  throat  requires  that  the 
coefficient  of  No  in  Equation  (2.21)  is  zero. 


(2.32) 


C  4  —  O 


_  Hr  Kr  t 

The  condition  that  and  are  continuous  at  the 

junction  of  Sections  A  and  B  yields  the  value 

i. 


(2.33) 


c*'- 


*(-  lrf:Y  v 


Hr  -%r  • 

The  condition  that  and  are  continuous  at  the 

junction  of  Sections  B  and  C  yields  the  values 
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where  (*  «~  Zra  )  Es/*  £j^(«.)  J^Ce.)  ■hJ_A.ce,>  J.^(€,)j 

_  a  —  -  *  X 


Clhla£l 


'«>*' 

where 


e0=  2.  E,^  X, 


By  the  use  of  the  previous  formulae,  the  expressions  a«/»-  is 
determined  completely.  The  quantities  K.  ,  R  ,  and  **  which 
determine  the  duct  geometry  are  determined  by  the  specific  problem 
under  consideration. 

Discussion  of  Results 

In  the  previous  formulae  the  quantities  ^  ,  and  **  are 

written  in  the  general  form.  This,  of  course,  enables  us  to  determine 
the  solutions  for  different  configurations. 

In  the  following  analysis  the  values  which  correspond  to  the  Wright— 
Patterson  Air  Force  Base  10— Foot  Y/ind  Tunnel  will  be  inserted.  Before 
this  is  done,  the  effect  of  the  choice  of  the  subsonic  junction  position 
X*  is  studied.  Most  investigations  of  a  supersonic  flow  are  made 
disregarding  any  influences  of  the  subsonic  portion  of  the  duct,  since 
there  is  no  theoretical  basis  for  its  consideration.  It  is  of  interest 
to  determine  some  numerical  values  of  the  subsonic  effects.  Therefore, 
a  flow  pattern  determination  will  be  made  of  a  configuration  which 
corresponds  to  the  supersonic  part  of  the  10- Foot  Wind  Tunnel,  and 
three  different  subsonic  junction  positions. 


Accordingly,  the  duct  shape  gives 

R  -  50  rQ 


oC  ■  .15  degrees 


Xq  ■  -0.10,  -0.20,  -0.U0 

The  values  of  <v)  are  the  Bessel  function  arguments  for  which 

I 


a^oand  are 


-  3.S32 


V"  7.015 
-10.1735 
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'Hie  -values  of  ,  £■  t-y.  ,  ^  -5^.  ,  and  E*.*  are  determined  by 
Equations  (2.12),  (2.13),  (2.17),  and  (2 .IS),  respectively.  The  junction 
constants  Cj^  ,  ,  and  C6«.  are  obtained  by  means  of  Equations 

(2.35),  (2.3 i|),  and  (2.35)*  respectively. 

Of  special  interest  is  the  velocity  distribution  along  the  axis  of 
the  duct.  The  x— velocity  terras  found  from  one— dimensional  theory  are 
given  ty  Equations  (2.3),  (2.6)  and  (2.9).  The  additional  x-velocity 
along  the  axis  caused  by  the  perturbation  potential  are  found  from 
Equation  (1.32).  Since  3%(«)=|,  one  obtains 

(2.56)  =  X  at; 

*  ' 

Figure  6  shows  the  *-<»,  values  through  Section  B  for  X0-  -,0i,-.o2  -.04,-00 
In  spite  of  appreciable  initial  flow  deviations  at  some  distance  upstream 
from  the  throat,  the  effects  are  hardly  noticeable  at  the  minimum  cross 
section  and  in  the  supersonic  region.  A  corresponding  behavior  is  found 
for  the  terms  ,  <***  ,  etc.  Accordingly,  even  extreme  changes  in 
the  subsonic  region  are  without  importance  for  the  supersonic  region. 

Figure  S  shows  the  velocity  distribution  along  the  axis  according 
to  the  one-dimensional  theory,  and  a  velocity  which  includes  the  correc¬ 
tions  given  ty  <*., ,  a*  ,  and  aj  .  Beyond  x  -  0.5  the  velocity  distri¬ 
bution  is  rather  smooth,  showing  the  present  configuration  to  be  satis¬ 
factory.  The  maximum  velocity  deviations  from  the  one-dimensional 
theory  have  a  magnitude  corresponding  to  a  Mach  number  of  about  .OQ[i. 

The  velocity  deviation  at  about  x  -  0.4  is  rather  large,  but  this  is  to 
be  expected  from  the  following  reasoning.  Point  M  (Figure  3)  is  con¬ 
sidered  to  be  the  junction  between  Sections  B  and  C.  The  effect  of  the 
change  of  shape  at  this  point  propagates  to  the  axis  along  the  Mach 
lines  MN.  Upstream  of  UN  the  flow,  especially  along  the  axis,  corresponds 
to  the  parabolic  surfaces.  Consequently,  any  change  in  pressure  distri¬ 
bution  caused  by  the  transition  is  felt  only  downstream  of  point  B. 

The  one-dimensional  theory  (Figure  g,  curve  <$>*„  )  shows  the  influence  of 
the  j -unctions  to  be  right  at  point  A.  The  superposition  of  the  <$* 
value  shows  the  occurrence  of  the  change  at  x  »  .38.  The  x-coordinate 
of  Point  N  in  Figure  3#  computed  by  an  exact  relationship  for  the  Mach 
angle,  is  x  •  .367.  Figure  S,  curve  (<9sK+$0  shows  what  might  be 
expected  from  a  more  general  reasoning.  The  functions  aQX'  , 

and  are  given  in  Figure  7. 

A  Remark  on  the  Application  to  a  Non-Symraetric  Configuration 

Although  no  examples  of  the  application  for  an  unsymmetric  flow  are 
carried  out,  an  idea  of  the  procedure  will  be  given.  The  shape  of  the 
duct  must  be  expressed  according  to  Equation  (1.6).  For  instance,  one 
must  carry  out  a  harmonic  analysis  of  the  function  r(x,»)  with  respect 
to  *0  for  various  values  of  x.  This  can  be  done  by  a  numerical  or 
analytical  evaluation  of  the  integrals  (1.6a)  and  (1.6b),  or  ty  means 
of  a  harmonic  analyzer  .  Then  the  values  of  are  found  from 

Equation  (1.9).  Following  this  preparation,  the  procedure  can  be  carried 
out  in  the  same  way  as  for  an  axially  symmetric  duct. 
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Solution  with  <Pe«  as  A  Function  of  % 

A  solution  of  the  differential  equation  is  determined  for  the 
coefficient  as  a  power  of  the  longitudinal  position.  The  solution 

and  its  derivation  will  be  found  in  the  Appendix  (Section  IV,  Part  F). 

Concluding  Remarks 

A  method  for  the  analysis  of  transonic  flow  through  a  nearly  circular 
duct  with  only  gradual  changes  in  cross-section  has  been  derived.  One  is 
finally  lead  to  the  solutions  of  ordinary  linear  differential  equations. 
(Equations  (1.33)  and  ( 1 . 3J4 ) ) •  These  solutions  can  be  carried  out 
numerically,  or  even  analytically  in  many  cases. 

Application  of  this  method  was  made  to  a  configuration  similar 
to  the  Wright-Patterson  Air  Force  Base  10-fbot  Wind  Tunnel.  A  velocity 
distribution  through  the  duct  was  determined  and  the  disturbances 
caused  by  the  variation  in  duct  shapes  were  computed.  The  velocity 
profile  for  the  configuration  is  satisfactory.  The  subsonic  influence 
in  this  case  is  small  enough  that  the  usual  consideration  of  no  subsonic 
disturbances  propagated  into  the  supersonic  region  is  justified. 
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SECTION  IV 


appenux  a 


Derivation  of  the  Potential  Equation 

It  follows  from  the  vorticity  theorems  that,  with  certain  condition 
for  the  initial  flow  which  are  fulfilled  in  most  practical  cases,  there 
exists  a  velocity  potential  <p  (X,  y,  a),  where  the  coordinates  are 
shown  in  Figure  1,  (That  is,  the  velocity  vector  can  be  expressed  as 
the  gradient  of  a  scalar  function  $  ) .  If  one  chooses  an  arbitraxy 
function,  <p  (X,  y.  A)  and  determines  the  velocity  components  accordingly; 

i«6i  j 

(5.1)  V*  =  <f>- 
(3*2)  —  <J>- 

(5.5)  Vw-A 

and  furthermore,  determines  the  pressure  from  Bermoulli's  equation,  then 
Euler's  equations  of  motion  are  automatically  fulfilled.  Therefore,  if 
a  velocity  potential  exists,  we  have  to  check  only  the  condition  of 
continuity.  The  continuity  equation  is  given  by 

^(?y  V,).^  *5(^  =  0 

If  one  inserts  Equations  (3.1),  (3.2),  and  (3.3)  into  Equation  (3.4) 


<J-5)  c %•)=© 

The  derivatives  of  p  are  found  from  Bermoulli's  equation  for  three- 
dimensional  flow 

(5.6)  A .  +  (<(>‘  *  0.'  -V  -  cov,sfoLvt 

where  i  denotes  the  enthalpy,  so  that 

'  —  . 

*  ~  *+ 1 

FYom  the  isentropic  flow  relations,  since  •*»  =  C y "T" for  perfect  gases. 

i-  -  (S-S?’1 

Ao  ) 

Consequently,  one  writes 


(5.7)  di  ~  c£  4®- 
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If  Equation  (3.6)  is  differentiated,  and  Equation  (3.7)  is  inserted 
the  resulting  expression  is 

c,‘  +■  <t>,  <t  <$>?  ♦ 

Or 

(3 JS)  =  - !> (<t>*  H*  +  <N  <1  <t>?  +  ** (-^f) ro 

Writing  Equation  (3.2)  in  detail,  one  obtains 


(3.2a) 


(3.2b) 


^?T=  4>y  <l>Sy 


.  ^U»  (^h»V 


(5.8c)  +  4>y  <t>^  +  ) 

Now,  if  the  partial  derivatives  of  ^  from  Equations  (3.2a),  (3*Sb), 
and  (3.2c)  are  inserted  in  Equation  (3«5)>  one  obtains  the  expression 

(3.9)  -  £;  O*  C p„  +  C(),  C^y  +  p^-)  +  ?  <t>« 

-ht®*  $ly  +  <Py  <h?  +  %T  + 

+  ^  <J>?y  -  fi  (Ix&ia.  * 


If  one  divides  Equation  (3*9)  by  ^  and  collects  the  terms,  one 
obtains  the  final  equation 

(3.io) 


This  is  the  complete  potential  equation  for  three-dimensional  flow. 
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APPENDIX  B 


Transonic  Law  of  Similarity 


The  general  three-dimensional  equation  for  a  steady  potential  flow 
is  given  by: 

*  $?/  0  -  )  +  <f><.  u 

-:>  <fcz<t>./<bgu,  <t>;  ' 

7* oj-  Y*-  a.1-  7 

The  velocity  of  sound  O-  is  obtained  from  Bernoulli's  equation. 


(3.12) 


.4!. ?  -  1=1 


Where  the  total  stream  velocity* 


N^J=C<l>^''-Ct)yX+  ^  * 


Let  us  consider  a  family  of  flow  patterns,  where  each  flow  pattern 
is  characterized  by  a  specific  value  of  a  parameter  2.  .  Let  the  velocity 
potential  be  given  by: 

0.13)  =  (x+ *<?(%,  y,  (5, y,u») 

where 

(3.i4) 

0.15)  'f-'j  2:  * 

This  hypothesis  is  justified,  by  the  fact  that  its  introduction  into 
the  differential  equation  for  $  ,  with  the  limiting  process  , 

results  into  a  non-trivial  equation  for  .  Indeed  substitution  of 

Equations  ("5.12),  (3.13),  (3«14)«  and  (3.15)  into  the  potential  equation 
(Equation  (3.11)1  gives 

(3.16)  a**  ($),-- (\-  <p _ /<£i.i+g$g)g*<py] 

*  Y v  <*-*-  t 

+  O*  Zr1  Cp^,  -  sT)  +  ^  Cp^  +  z0 


to.  i‘  cPr  . 

77  cx.  *>  /  1  Y  y  x 

The  value  of  the  sonic  velocity  as  given  in  Equation  (3.12)  dan 
be  written  in  terms  of  Or  and  the  velocity  potential,  and  is  finally 
expressed  as 

(3.17)  OL**  a* 
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The  introduction  of  Equation  (3*17)  into  Equation  (3.16)  leads  to 
(3  *1&  J  ^ 

Li  -(*->)£  ‘P*  D  -(*->)£<¥?] 

+  £l  D-Cn-i)a^]  =  o 

Since  2:  is  small,  only  the  terms  of  the  lowest  power  of  £  prevail. 
The  equation  is  divided  by  ,  and  one  obtains 

(3.i9)  -(*+■'>  cp«  cs?7  +  q>Ty  *  ^  =  o 

This  equation  shall  be  referred  to  as  the  simplified  transonic  potential 
equation. 

Let  us  consider  Equations  (3«13)>  (3.1 h)>  and  ( 3 » 15 )  for  small 
values  of 

C3.20)  $  -  0.*i (  *  i-  2:  Cp  ( K, Xft  U>1) 

The  velocity  components  given  by 

(3.21)  <$tr  (  x,  y,  **0 

(3.22)  -  gJ*2:  X  <%  Cx}y,^ 

(3.23)  2^ 

7  7 

Let  vis  consider  flow  patterns  which  belong  to  different  values  of 
the  parameter  z  .  In  such  flow  patterns  the  points  having  the  same  x 
and  y  values  are  called  corresponding  points.  From  Equations  (3*14) 
and  (3*15)  it  is  seen  that  at  corresponding  points,  the  values  of  H 
are  the  same  and  the  y  values  vary  as  h’*.  At  corresponding  points 
the  deviation  of  the  x  velocity  is  proportional  to  Z  ,  the  y  velocity 
and  the  <*>  velocity  are  both  proportional  to  .  The  significance  of 
*  is  found  from  the  conditions  at  infinity.  At  infinity  the 
deviation  of  the  X  velocity  from  the  sonic  velocity  is  proportional  to 
the  difference  of  the  free  stream  Mach  number  from  one.  Therefore, 

2  is  also  proportional  to  the  difference  between  the  stream  llach 
number  and  one. 

For  Z:  -  0  the  streamlines  have  a  constant  y.  The  deviations  of 
the  streamlines  from  a  line  (  7  “  constant)  are  found  by  integrating  the 
slope  of  the  streamlines.  The  component  of  the  velocity  deviations  oft**, 
as  well  as  the  deviation  of  7,  behaves  proportional  to  •  Conse¬ 

quently,  the  deviation  of  the^streamlines  from  the  lines  y  ■  constant, 

<*»  --  constant,  behave  like  a  t  .  To  summarise  the  transonic  law  of 
similarity: 
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1.  The  x  coordinates  of  corresponding  points  are  the  same 

2.  The  y  coordinates  of  corresponding  points  are  proportional 
to  2"t 

3.  The  difference  between  the  x  velocity  and  the  sonic  velocity 
is  proportional  to  2. 

The  difference  between  the  pressure  and  the  sonic  pressure  is 
proportional  to  2 

5.  The  deviation  of  the  streamlines  from  lines  of  constant  y  and 
u»  is  proportional  to  •£  T 


APPENDIX  C 


Boundary  Condition 

Along  a  booty-  at  rest  the  velocity  vector  is  parallel  to  the  surface; 
that  is,  the  component  of  flow  normal  to  the  surface  is  zero.  To  express 
the  relation,  let  the  unit  vectors  in  the  direction  of  increasing  x,  y 
and  <*»  coordinates  be  Ux,  Uy,  and  U  w  ,  respectively.  The  total 
velocity  vector  then  is  expressed  by 

(3.2ii)  <£[C\+  <P0  tA*  +  Cpy  Uy  +  Um] 

The  vector  component  normal  to  the  duct  surface  is  given  by 

(3.25)  "  W-x  ~  +Uy  =  0 

The  scalar  product  of  these  two  terms  is  the  component  of  the  velo¬ 
city  vector  normal  to  the  surface.  According  to  the  boundary  condition, 
this  component  is  zero.  Then  one  obtains 


(3.26)  CP*>  +  of  ^  =  0 

The  assumption  that  $  is  a  small  quantity  implies  that  the 
slope  of  the  duct  and  surface  is  gradual,  and  consequently  ^L.  and 
are  small  quantities. 

Therefore,  the  expressions  —  and  ll-  are  of 

*  in  ' 

second  order,  and  can  be  neglected  in  Equation  (3.26).  This  results  in 


(3.27)  - 


Vv 
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APPENDIX  D 


Derivation  of  the  Basic  flow  Velocity 


The  basic  flow  velocity  is  the  velocity  obtained  from  one- 

dimensional  flow  theory,  and  is  expressed  as  a  function  of  the  duct  radius 
at  the  cross  section  under  consideration* 

The  continuity  equation  is  written 

(3.28)  r*  tt  p  V  =  r*  tt  p  <*• 

where  V  is  the  absolute  value  of  the  velocity.  Since  is  a  function 

of  V.  it  can  be  developed  with  respect  to  (V  -  a#).  une  then  obtains 

*■  o*  ,  v  (V-o*)1 

(3.29)  =  P  - £ 

If  one  introduces  the  relation 

(3.30)  A  r  -  Y*—  Y^0 

the  continuity  equation  is  rewritten  as  - 

r^of-  ]-r(n+w>,*e  °*  Tr'r-1 

If  r  is  small,  the  higher  order  terms  can  be  neglected  so  that 

(  *  \  -  r  Ar  ~~|  i 

0.51)  CV-6.  )  -\_V5T  -y*J 

Since  ,  the  basic  flow  velocity  is  expressed  by 

v  or 

(3-52)  n-]* 

APPENDIX  E 

Analytical  Representation  of  the  Constants  T.^,  ,TVm.  ,  and 


An  analytical  representation  is  sought  for  the  integrals  that  occur 
in  the  expressions  and  Th®  first  term  X*^is 

(3.33)  To^=  f  y  ^o^C7)<4y-  f'y  ^ 

*  J«  Jo 

Integrating^  _  y*| 

(5.3U)  1.^-  £  D.'Civ^VS.  (v«>*  v)  Yur. 


3ibid.,  p.  Hi5 
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Then**  since 

T-,  y)  =  7,  y) 

(3-55)  =  t  l^C  V v)  +  V  ^  y)J 

The  values  for  the  arguments  of  the  Bessel  function  7.  Cv  y)  were 
chosen  so  that  *  0.  However,  it  is  also  true  that  ^T„ (~y) 

so  that  Equation  ($.35)  simplifies  to  the  analytical  expression  ^**1 

The  term  Koj*.  is  defined  by 

(3.37)  K.r  =  J^'  ys  2-.^Cy)  4y  =  y s  T. C%  y)  iy 

If  one  denotes  the  argument  by  ^  =  )  Equation  (3  *37)  is  rewrit¬ 
ten  as  ' 


or  simply 


(3.3S)  K.  =  *  U.  "50  A  vx. 

/or 

An  integration^  gives  the  relationship 
(3*39)  K^s 

or  ° 

«•'“»  K„^=  (.)u  u‘ J.Cu)  -  4U  T,(^)].V 

If  the  Equation  (3.A|0)  is  rewritten  in  terms  of  y,  and  the  limits 
of  y  **  0  and  y  ■  1  are  inserted,  the  new  relation  becomes 


(3J4l) 

Since 

(3J42) 


K“r-  t;  t.2-  ^  7.  ( i  y)  +  Cv  -  *  i.r)  T,  (54  r)] 

— ^  ^  =  =  O  the  \(*jk  is  reduced  to 


K.^=  "S.Ovy) 


**  Ibid.,  p.  li*5 
5  Ibid.,  p.  145 
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The  term  1®  defined  as 

(3.U3)  ^CyMy*  £  v  V(%.  v)Ay 

It  is  seen  that  the  integration6 7  leads  to  the  expression 

iiM  xvr=  ^  [.V  (%r  y)-1vH  <%.  y)1*«  (Vry)X 

Also7 

(3.U5) 

(346)  -  Zr+1  (7)  -  ^r1'1  --Sf-ifCy) 

The  substitution  of  Equations  (3.45)  e-nd  (3 .46)  into  Equation  (3.44) 
gives  the  new  relation 

(347)  £  tVCvY)+  V 


where  Tv#  denotes  the  firwt  derivative  of  with  respect 

to  the  argument. 

Since  the  arguments  are  chosen  such  that  3" v  C  "*)  vj.  0  and 

the  integration  limits  are  y  ■  0  and  y  »  1,  the  value  of  Xv^.  is 

(3.66) 

The  term  K  vj*.  is  written 

(349)  K.r= JV+’  s.V'*')  Ay =  [' 

Substituting  the  following  differential  equation 


(3.50)  ^  Cy  ^v^v^y\1  ~  y  ]  Y^ 


into  Equation  (3*49) ,  one  obtains 

(3.51)  Kv^=  3^  {'S.Y''  *7  fr'5*  ^3rViJAy 


6  Ibid.,  p.  146 

7  Ibid.,  p.  145 
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Integrating  by  parts,  one  can  rewrite  Equation  (3*51) 


0.52)  K,^=  ~  [-/yT^v  r)J/^yvjJ  (%) 4 y + v‘jfv ,w  XcVjMy] 

Further  integration  and  the  substitution  of  the  limits  of  y  -  0 
and  y  ■  1,  leads  to  the  relation 

(j.53)  K,^=  ^Crlv  (^.ry)+  vTv^v^y)] 

Since  y^sothe  analytical  representation  of  VC -y  j*.  ig 

APPENDIX  F 

Solution  of  Equation  (1.33)  if  a  *  is  Given  by  X  * 

The  solution  of  the  differential  equation  (1,33)  is  determined 
for  a  more  general  case.  Let  the  coefficients  of  the  differential  equation 
be  such  that  it  is  a  power  t  of  the  longitudinal  position. 

(3*55) 

This,  when  inserted  into  the  differential  aquation  gives 

(3.56)  cl  y  a*  -v  — *  (  a.*  -  o 

In  order  to  remove  the  inconvenient  power  X  in  the  coefficient 
°f  i  1st  us  introduce  a  new  independent  variable  u.  by  the 

equation  u  -  X  where  k  is  a  suitable  constant  which  is  chosen  later. 

Then  one  obtains 


(3.57) 


(3.5S) 


ii-r  =.  |(  v1"'  4  o£ 


du. 

■  %  X 


If  one  writes  X-  and  inserts  Equations  (3. 57)  and  (3.5S) 

into  Equation  (3.56),  the  following  equation  will  result 


(3.59) 


»l  *  1  #  x  * 

~7Z?r  *  Tv  -5^-  +  v. 
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where 


(3.60)  r  = 

If  ^  s  |  -  2L.  t  Equation  (3.59)  simplifies  to 

(3.61)  +  i-  £  *  tK  q.%  =  o 

du1  ^  du.  " 

Equation  (3.61)  has  the  following  form 

(5.62)  Cl."  +  (  O-'  +  (V  +  ^  °  . 

s 

the  solution  of  which  is  given  by 

(3.63)  C^=US£p(^  ia) 

where  By  is  a  Bessel  function  of  p  order. 

let  us  set  the  following  relations: 

(5.6U)  *  = 

(5.65)  p  -  *  = 

(3.66)  ^  ~  (jT5) 

The  solution  of  the  differential  equation  (3*61)  involving  coeffi¬ 
cients  that  are  powers  of  x  is  expressed  as 

(5.67)  ~  ^*£(.v)Zul 

where 

(3.66)  S  =  £§ 


8  Ibid.,  p.  346 
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G  I  —  CYLINDRICAL 
COORDINATE  SYSTEM. 


FIG. 2-  SPACE  ELEMENT  IN 
CYLINDRICAL  COORDINATE  SYSTEM. 


FIG.  3“  PROPAGATION  OF  DISTURBANCES  FROM  SUPERSONIC  JUNCTION 


FIG-4—  AXIALLY—  SYMMETRIC  DUCT  CONTOUR 


AMC  ffilOHT  KEELD 

AIRCRAFT  LABO.uTORY  wIND  TUNNEL  BRANCH 
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Figure  6.  -  Mach  Number  Deviation  (From  One-Dimensional  Theo; 
for  Different  Subsonic  Ducts. 


£MC  YffilGHT  FIELD 

AIRCRAFT  LABORATORY  WIND  TUNNEL  BRANCH 


Figure  S.  -  Additional  Mach  Number  Along  the  X-Axis 


